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, High rank solutions to the 2D Toda Lattice System are constructed simultaneously 

with the effective calculation of coefficients of the higher rank commuting ordinary 
difference operators. Our technic is based on the study of discrete dynamics of Tyurin 
i parameters characterizing the stable holomorphic vector bundles over the algebraic 

0^ ' curves (Riemann Surfaces). 

o 

^ : 1 Introduction, 
o 

. The aim of this work is to consider a circle of problems closely related to the construction 
of algebro-geometric high rank solutions of the two-dimensional (2D) Toda lattice 



^ : df^ipn = e'^"-'^"-^ - e'^"+i-^". (1.1) 

High rank solutions of the Kadomtsev-Petviashvili (KP) equation 

X ' SUyy = {AUt - QUUx + Uxxx)^ (l-2) 



were constructed in the earlier work of the authors Their construction, as well as the 
notion of rank of the solutions was based on the theory of commuting ordinary differential 
operators |2]. In modern mathematical physics this theory arisen as an algebraic by-product 
of the integration theory of soliton equations and the spectral theory of periodic finite-gap 
linear differential operators 1311113111111 initiated in 

The KP equation, as well as any other soliton equation, is a part of the corresponding 
hierarchy of commuting nonlinear partial differential equations. A solution u = u(t) of the 
KP hierarchy is a function of an infinite set of variables t = {ti = x,t2 = y^t^ = t,t4, . . .). 
Beginning with the work [Zj, the algebro-geometric solutions of integrable (l+l)-systems of 
the KdV type are singled out by the constraint that they are stationary solutions for one 
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of the flows of corresponding hierarchy. The algebro-geometric solutions of the KP equation 
were found in the work They are singled out by the constraint that they are stationary 
for two flows of the hierarchy dt„u = dt^u = 0. This constraint is equivalent to the existence 
of a pair of commuting ordinary differential operators 

n m 

[Ln, Lm] = 0, Ln = Y.'^i^x^ ^rn = J^'^J^l ' i^'^) 
i=0 j=0 

that commute also with the auxiliary Lax operators dy — L2, dt — for the KP equation 
found in 013, 

L2 = dl + u, Ls = dl + ^ud, + w. (1.4) 

Therefore, pairs of commuting operators define invariant subspaces for the KP hierarchy. 
By definition, rank r of a pair of commuting operators is the number of linear independent 
common eigenfunctions, i.e. the number of solutions of the equations 

Lnilj' = E^\ L^^P' = w^\ ^ = l,...,r, (1.5) 

for {E, w) such that a solution exists. The rank of the pair of commuting differential operators 
is a common divisor of orders of the operators. 

The classification problem of commuting ordinary differential operators with scalar coeffi- 
cients was considered as a pure algebraic problem by Burchnal and Chaundy in the 20-s 
[Tfl| ITT] . They solved the problem for the case of operators of co-prime orders (where the 
rank always equals 1) although no explicit expressions for the coefficients of the operators 
were found. Effective classification of commuting rank 1 operators were obtained in [Ij. 
Burchnal and Chaundy emphasized (see [H]), that the problem for the case of rank r > 1 is 
very complicated. 

The first steps were made in [121 HHI- A method of effective classification of commuting 
generic differential operators of rank r > 1 was developed by the authors in |U[2]. Commuting 
pairs of differential operators rank r depend on (r — 1) arbitrary functions of one variable, 
smooth algebraic curve T with one puncture P and a set of the Tyurin parameters for a 
framed stable holomorphic vector bundle over T. The corresponding constructions are called 
one-point constructions. 

For the 2D Toda lattice equations the constraints that single out algebro-geometric 
solutions can be described in the same way as in the KP theory. The condition that these 
solutions are stationary for two flows of the hierarchy is equivalent to the existence of a pair 
of commuting ^, //-dependent difference operators 

JV+ A/+ 

L= <ri)T\ A= <ri)T\ (1.6) 

i=-N- i=-M- 

that commute also with the Lax operators 

U = d^-T-win), U = d,, - c{n)T-\ (1.7) 
w{n) = ipn^, c(n) = e^"-^"+S (1.8) 
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for the 2D Toda equations. Here and below T denotes the shift operator, Ti/n = Un+i- 

The well-known theory of commuting rank r = 1 scalar difference operators was made 
of only two-point constructions jHJUn]. A ring of such operators is isomorphic to the ring 
A(r, P^) of meromorphic functions on an algebraic curve F with the poles at two punctures 
P^. Remarkable, almost graded structure of these rings, their connections with the Fourier- 
Laurent type basis on Riemann surfaces, was developed by the authors and used for the 
string theory [T6] . 

The classification problem of commuting difference operators in its full generality has not 
been solved yet. Analysis of the problem presented in section 2 of this article has allowed 
the authors to obtain some new important results. These results were announced in [^1^ . 
It turns out that commuting rank r = 1 difference operators can be obtained with the 
help of multi-point constructions with any number of infinite points including one point if 
rank is even (which is most astonishing). Multi-point constructions are well known in the 
theory of commuting differential operators with matrix coefficients It seems reasonable 
to consider commuting difference operators constructed with the help of only one- and two- 
point constructions, as a difference analog of commuting differential operators with scalar 
coefficients. For the authors the supporting evidence of this suggestion is that only in these 
two cases the rings of commuting operators are invariant with respect to 2D Toda hierarchy. 

It is necessary to emphasize that one- and two-point constructions differ drastically. 
There are no functional parameters in the two-point case. Coefficients of the corresponding 
operators can be written in terms of the Riemann theta- functions. The one-point construction 
of commuting operators of rank / > 1 contains functional parameters. The coefficients of 
commuting operators depend on I arbitrary functions of discrete variable, smooth algebraic 
curve F with puncture P and a set of Tyurin parameters. As in continuous case [2], the 
reconstruction problem of the eigenfunctions of the corresponding operators is equivalent to 
the Riemann problem on F, and in general can not be solved explicitly. At the same time, 
as it was found in |T], in some cases the reconstruction problem for the coefficients of the 
commuting operators can be solved explicitly. The solution is based on equations for the 
Tyurin parameters describing their dependence on normalization point. Note, that in the 
recent work of one of the authors connections of continuous deformations of the Tyurin 
parameters, describing high rank solutions of the KP equation, and the Hitchin systems 
were established [IB]. Using a discrete dynamics of the Tyurin parameters, the authors 
found explicit form of rank 2 commuting operators corresponding to an elliptic curve with 
one puncture. The ring of such operators is generated by a pair of operators L and A of 
orders 4 and 6, respectively: 

L = ^ u,{n)T\ A=j2 Vi{n)T\ 

i=-2 i=-3 

The coefficients of the operators depend on two arbitrary functions 7„ and s„. In the simplest 
case an explicit form of L is given by the formulae: 

where L2 is the difference Schrodinger operator 

L2 = T + Vn + CnT~'^ 
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with the coefficients 

2vn+i = s„F(7„+i,7„) - s„+iF(7„,7„+i). 

where 

F{u, v) = Ciu + v)- Ciu -v)- 2Civ). 

Here and below p{u) = p{u\uj,u!'), ({u) = ({u\uj,u') are classical Weiershtrass functions. 

These formulae are discrete analogs of the formulae for differential operator L4 of the 
order 4 and rank 2, obtained by the authors in [12]. In [20] an explicit expression of the 
functional parameter in the formulae for L4 corresponding to the Dixmier operator |12| was 
found. Coefficients of this operator are polynomial functions of the independent variable x. 
It would be interesting to find a discrete analog of the Dixmier operator. 

The construction of high rank algebro-geometric solutions of the 2D Toda lattice was 
briefly described in [2I1I221- We discuss it in detail in Section 3. In is based on the construction 
of the vector Baker- Akhiezer function as a deformation of the eigenfunctions of commuting 
differential operators. These deformations are defined by the form of the Baker-Akhiezer 
functions in the neighborhood of punctures. The form is defined with the help of some germ 
matrix functions. Once again it is necessary to stress the difference between one and two 
puncture cases. 

In the two-point case (in which there are no functional parameters in the construction 
of commuting differential operators), the Baker-Akhiezer function is defined by two auxil- 
iary germ functions z). Each of them is defined by an ordinary differential 
equation containing arbitrary functions of one of the continuous variables ^ ov rj, respective- 
ly. In the one-puncture case, functional parameters are present in the construction of the 
corresponding commuting operators, but there is no ambiguity in the definition of the germ 
function \l/o(n, ^, 77, 2;). 

In both the cases the multi-parametric Baker-Akhiezer functions ipn = (ipn) satisfy the 
equations 

LiVn = 0, L2^„ = 0, (1.9) 

where Li, L2 are operators of the form ()1.7jl . Compatibility of (ll.9jl is equivalent to (ll.ljl . 
Equations HIM imply also that for each n the Baker-Akhiezer function, as a function of the 
variables (^,^7), satisfies the linear Schrodinger equation: 

+ ^«(^' V)d^ + «"(^' V)) = 0. (1.10) 

Therefore, our constructions provide at the same time the construction of two-dimensional 
Schrodinger operators in the magnetic field that are integrable on one energy level. 

Inverse spectral problem on one energy level for the Schrodinger operators in a magnetic 
field was introduced and solved in [231 the rank 1 case. The construction oi2D Schrodinger 
operators of arbitrary rank integrable on one energy level was proposed in ^H]- All these 
constructions are two-point constructions. The possibility of the one-point construction of 
integrable Schrodinger operators has never been discussed until now . 
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We would like to stress that the problem of construction of solutions of the 2D Toda 
lattice equations is equivalent to the problem of construction of integrable Schrodinger op- 
erators in a magnetic field. In the modern theory of integrable systems equations (jl-Hl and 
its Lax representation were obtained in the framework of the Zakharov-Shabat scheme as 
two dimensional analog of the one-dimensional Toda lattice ([21]). Later it became clear 
that these equations had been known in the classical differential geometry in an equivalent 
form of chains of the Laplace transformations for the two-dimensional Schrodinger operator 

(I2ni)- 

Let us consider the two-dimensional Schrodinger operator L in a magnetic field 

2L = {d, + B){d, + A) + 2V, (1.11) 

where = — idy, 0^ = 0^ + idy. Usually the functions H = {B^ — Ag)/2 n U = —H — V 
are called magnetic field and potential, respectively. Below, for brevity, the function V 
would be called the potential of L. The operator L is defined up to the gauge transformation 
L e^^ Le^ . The only invariants of L are the potential V and the magnetic field H. Laplace 
transformation is a two-dimensional analog of the Baclund-Darboux transformation. It is 
defined as follows 

{d-, + B){d, + A) + 2V ^ V{d, + A)V-\d, + B) + 2V. (1.12) 

The potential V and the magnetic field H of the transformed operator are given by the 
formulae 

V = V + H, 2H = 2H + A\nV. (1.13) 

If a function ip satisfies the equation Lip = 0, then the function -ip = {dz + A) ip is a solution 
of the equation L ijj = 0. Let us consider a chain the Laplace transformations 

Vn+l = Vn + Hn+l, 2Hn+l = 2Hn + AlnVn. (1.14) 

After the change of variables Vk = In (</)„ — fk-i) the chain of the Laplace transformations 
becomes equivalent to the the 2D Toda lattice equations (ILljl . 

2 Commuting difference operators 
2.1 Formulation of the problem. 

The main goal of this section is to construct an effective classification of commuting 

[L,A]=0 (2.1) 

scalar difference operators of the form 

N+ M+ 

L= Ui{n)T\ A= N± = r±n±, M± = r±m±, (2.2) 

i=-N- i=-M- 
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where r± are maximum common divisors of maximum and minimum orders of the operators, 
respectively, i.e. the non-zero numbers n±,m± are co-prime: 

(n+,m_,.) = {n_,m_) = 1. (2.3) 

Let us assume, that the leading coefficients of these operators are different from zero. For 
brevity, they are denoted by 

M±(n) = u±N±{n) ^ 0, v^{n) = v±M±{n) ^ 0. (2.4) 

Equation ()2.H) is invariant under the gauge transformation 



L^A ^L = gLg-\A = gAg'\ L = J2 g{n + i)g-\n)ui{n)T\ (2.5) 

i=-Af_ 

where g{n) 7^ is an arbitrary non-vanishing function of the discrete variable n. This 
transformation can be used for the following normalization of L 

M+(n) = 1. (2.6) 

The normalization ()2.6p which will be always assumed below, is invariant under the gauge 
transformations corresponding to the functions g{n) such that g{n + iV^) = g{n). Equation 
(12.111 implies that the leading coefficient of the operator A satisfies the equation v^{n+Nj^) = 
v^{n). Therefore, the gauge transformations can be used in order to fix the normalization 

v'^in + r+) = v^{n) = (2.7) 

in addition to ()2.6jl . Here n is a residue of n mod r_|_, i.e. n ^ n = n(mod r_|_), < n < r_|_ — 1. 
The normalization ()2.6l2.7jl is invariant under the gauge transformations corresponding to 
the functions g{n) such that g{n + r+) = g{n). 

If the normalization (j2.6l2.7jl is fixed, then for the coefficients on the other edge have the 
form 

u~{n) = h'^{n - N_)h{n), v~{n) = v^h~^{n - M_^)h{n). (2.8) 
Here n = n(mod r ). 

In the case r+ = r_ = r > 1 the commutativity equations have additional symmetry. 
Namely, let ?/„ be a function of the discrete variable n, then for any < i < r we define a 
new function Yn 

jo, n^^imodr) 
[ i/n, n = rn + I 

Under that correspondence a difference operator of order defines an operator of order 
rN. The direct sum of r pairs of commuting operators [Li, Ai] = considered as opera- 
tors acting on defined above sublattices, satisfies a non-degeneracy conditions (12. 4|) and the 
commutativity equation ((231) • 

Below we will consider only commuting pairs of irreducible operators, i.e. operators that 
have no block-diagonal or jordan-cell form under decomposition of the n-lattice into a sum 
of sublattices kr + i, 0<i<r — 1. A pair of commuting operators is irreducible if, for 
example, for each i = 1, . . . ,r — 1 there is tiq so that 

Vrm+-i{no) 7^ 0, V_rm-+iino) ^ 0. (2.10) 
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2.2 Spectral curve. Formal infinities. 

In the discrete case an affine spectral curve is defined identically to the continuous case. Let 
us consider a linear space C{E) of solutions of the equation Ly = Ey, i.e. 

N+-1 

yN+ + Ui{n)yn+i = Eyn. (2.11) 

i=-Af_ 

Dimension of this space is equal to the order of L, dmiC{E) = + iV_. Restriction of A 
onto C{E), 

A{E)=A\ciE), (2.12) 

is a finite-dimensional linear operator. The spectral curve, which parameterizes common 
eigenfunctions of L and A, is defined by the characteristic equation 

R{w,E) = det{w - A{E)) = 0. (2.13) 

Let c* be a basis of solutions of (|2.1H1 . defined by the initial conditions 

ci = Sl i,n = -N^,...,N+-l. (2.14) 

Matrix entries of A{E) in that basis are polynomial functions of the variable E. Therefore, 
R{w, E) is a polynomial in both the variables w and E. 

Compactification of the spectral curve. In the discrete case a construction of common 
eigenfunctions of commuting operators at infinity — *• cxd is on the whole parallel to the 
continuous case, but at the same time contains some new important features. 

Let L_)_ be the linear space of solutions of the equation 

L<l> = 2-^+$, $ = {$„} (2.15) 

of the form 

^n{z) = z~'' (Y.Un)z'^ ■ (2.16) 

It is assumed that ^o("') is not equal to zero for at least one value of n. 

Lemma 2.1 The space considered as a linear space over the field fc+ of Laurent series 
in the variable z, is of dimension A^+. It is spanned by the solutions $*, i = 0, . . . , — 1, 
uniquely defined by the conditions 

$^(z)=z-"5„„ n,^ = 0,...,iV+-l. (2.17) 

For the proof of this Lemma it is enough to note that a substitution of the series (j2.16ll into 
(I2.15jl gives a recurrent system of equations for the coefficients of (I2.16jl 

Un + N+) = Un), 

Un + N+) = ^i{n)-Un+N+-i{nMn + N+-l), ... (2.18) 
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The system implies that $ is uniquely defined by the initial data Csi^), i = 0, . . . , A^_|_ — 1. 
The space L_|_ is invariant with respect to the operator A. Therefore, 

A^+-l oo 

A^niz) = E M^Mi^)^ M^) = E 4^^- (2-19) 

j=0 s=-M+ 

The leading coefficient is a matrix 

A^j^^^ = V^6j^i+M+{mod N+) (2.20) 

The change z ez, where = 1, defines an automorphism of L+iz). That implies that 
coefficients of the characteristic polynomial 

det{w - Aij{z)) =w^+ + J2 cii{E)w\ E = z-^+ , (2.21) 

1=0 

are series in the variable E = z~^+. From that it follows that the eigenvalues are Laurent 
polynomials in the variable z (but not in fractional powers z^^'', k > 1). 

Let us assume that 

vtl^vl, i^j, (2.22) 

where are defined in (|2.7|1 . In that case the matrix Aij{z) has a unique eigenvector 
^Wf^), i = r, - 1, of the form (ET 



^^\z) = Z-- l^6,,n + T.UnVj . (2.23) 

The leading term of the expansion is an eigenvector of the matrix (I2.2()L corresponding to its 
eigenvalue . Under the change z — > e'^+z, the leading term gets the factor Therefore, 
the vector z*\Ef^*''(z) has an expansion in powers 

^0 = = (2.24) 
That implies for 2 > and j = 0, . . . , r_|_ — 1, the following equalities: 

^£++,(^) = z-'z,\Oizo)), j<z, 

^il+,(^) = z-^z,\0{l)), j>t. (2.25) 

From the irreducibility condition ()2.1()jl . in which without loss of generality we put uq = 0, 
it follows that for i > the zero component of the eigenvector has the form \l/o^ = 0(2;'""*). 
Therefore, the normalized eigenvectors 

^1^^ = ^, # = 1, (2.26) 

n\z) 
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for i > have the form 



V^£++,(^o) = Oiz,"), j<i. (2.27) 

For i = we have: 

V'£U(^o)=^o~'(l + 0(^o)). (2.28) 

Hence, in the case when the conditions ()2.22|) are satisfied, we have constructed a set of 
formal eigenvectors of the operator A{E) in the neighborhood of the infinity E = 00. Below, 
we construct an additional set of formal eigenfunctions under the assumption 



where f , are defined in (E 



vr^^i, ^^J, (2.29) 



Let us consider a linear space L_ over the field of Laurent series k_ in the variable z_, 
spanned by the solutions of the equation 

L$ = z-^~ $ (2.30) 

of the form 

Ki^) = ^-l^f:Cin)z!_y (2.31) 

It is assumed that ^o{n) is not equal to zero for at least one n. 

As before, we obtain that the characteristic polynomial of the operator A~{z_), induced 
by an action of A on L_, has the form 

r_ n_ 

det{w-A-{z.)) = l[l[{w-vr{zk,-)), z^,- = e^zL, (2.32) 

i=0 k=0 

where e"" = 1, and the series = v^{zq^_) has the form: 

viizo,^) = i;r^o:""(l + 0{zo,-)) . (2.33) 

The normalized eigenvector of A^{z_), corresponding to the eigenvalue w = v~{zq^_) for all 
z = 0, . . . , r — 1, has the form: 

^S.+jM = o(4_), t>j, 

V'il+.l^-) = 0(4+1), ^<J. (2.34) 

In the direct sum of the spaces L+ and L_ it is possible to choose a basis c^, normalized by 
the conditions 

4 = (5i,„, z,n = N_,...,N+-l. (2.35) 
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In this basis the operator A has the same matrix entries as in the construction above of the 
affine spectral curve. Therefore, the characteristic equation ()2.13|) coincides with the product 
of the characteristic equations for A{z) and A~{z_), i.e. 



det{w-A{E)) 



r+ — 1 

n n {w-vi{zk)) 

i=0 k=0 



T- n— — l 

n n {^-K{zk,-\ 



j=0 k=0 



(2.36) 



where the products in the first and the second sets of factors are taken over all roots of E ^ 
of degree n+ and n_, respectively: 



E 



(2.37) 



The expansion ()2.36^) gives a complete information on compactification of the spectral curve 
when the operator A{E) has N distinct eigenvalues for almost all E. It is the case of com- 
muting operators of rank 1, which is considered in the next section. 



2.3 Commuting operators of rank 1. 

The conditions (I2.22jl are ()2.29jl imply that the eigenvalues of A{E)^ corresponding to factors 
of the first and in the second products in the formula (I2.36jl are distinct. In addition, let us 
require that the factors in the two products do not coincide with each other. For that it is 
sufficient to assume one of the following conditions 

(i) m+n_ 7^ m_n+, {ii) vf ^ vj . (2.38) 

In this case equation (I2.36jl implies that the affine spectral curve defined by equation (I2.13jl , 
is compactified by / = (r+ + r_) points P^, i± = 0, . . . , r± — 1. In the neighborhood 
of the infinity, and therefore for almost all E the spectral curve T has N = + N_ 
sheets. Therefore, to every point of the spectral curve there correspond the unique up to 
multiplication common eigenfunction ipn of the operators L and A. Let us present the main 
theorem. 



Theorem 2.1 Let irreducible commuting operators satisfy the conditions ^2.2i^[^7^\2.S^} . 
Then: 

(1) the spectral curve V, given by the characteristic equation ^2.1,^) is compactified at 
the infinity by I = (r_|_ + r_) points P^, in the neighborhoods of which local coordinates are: 

Zk,± = E-^l^^- 

(2) the common eigenfunction of the commuting operators 

LMQ) = EMQ), AMQ) = wMQ), Q = {w,E)er, 

normalized by the condition ipo = 1 is a meromorphic function on V , with the divisor of 
poles D = {7s} outside the punctures P^ does not depend on n. In the neighborhoods 
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of the punctures ipn has the form ^2.2'T^2. 2^2. 34\) , ie. if n is represented in the form 
n = kr+ + j+ = k'r_ + j-, < j± < r±, then 

(2a) ipn has the poles of order k at the punctures Pfl, Pj^+i, • • • , P^_^_-i and the poles of 
order k + 1 at the punctures Pi, . . . , Pj^; 

(2b) ipn has zeros of order k' at the punctures Pj~_, . . . , P,r _i and zeros of order k' + 1 
at the punctures Pq, . . . , Pji-i 

(3) in the general position, when the spectral curve is smooth and irreducible, the number 
of poles 7s (counting with their multiplicities) ofipniQ) outside the punctures is equal to the 
genus g of the curve T. 

The theorem defines a map which for a generic pair of commuting operators, satisfying the 
conditions of the theorem, assigns a smooth algebraic curve T with / punctures P^ and a 
divisor of degree g: 

[L,A] = 0^{T,Pt^,D = {^s}}, 0<t±<r±, s = l,...,g. (2.39) 
Let us show that these algebro-geometric data define uniquely the commuting operators. 

Let r be a smooth genus g algebraic curve with / = (r+ + r_) punctures P^. From the 
Riemann-Roch theorem it follows that for any non-special divisor D = (71, ... , 7^) there is a 
function ipn{Q) unique up to a factor, such that its divisor of the poles outside the punctures 
is not greater than and such that at the punctures P^ it has poles and zeros of the orders 
defined in (2a, 26). 

Indeed, the conditions (2a, 26) mean that ipn{Q) belongs to a space of meromorphic 
functions L(D„), associated with the divisor Dn 

r+ — 1 jj^ r_ — 1 

i+=0 1=1 i_=0 i-=0 

where n = kr+ + j+ = k'r_ + j_. This divisor is of degree g, therefore, according to the 
Riemann-Roch theorem the space L(D„) is one-dimensional. 

Let us denote a ring of the meromorphic functions on T with poles at P^ by .4,(r, P^). 

Theorem 2.2 Let ip{Q) = {ipn{Q)} be a sequence of functions, corresponding to algebro- 
geometric data ^2.3fil\) . Then, for each function f G A{T, P^) there is a unique difference 
operator Lf (with coefficients independent ofQ), such that 

Lfi;{Q) = fiQ)i;iQ). 

If the function f{Q) has poles of order and n_ at the punctures P^ respectively, then 
the operator Lf has the form ^2.2^) . 

The proof follows immediately from the Riemann-Roch theorem. If the function / has poles 
at Pj^ of order n±, then the function f{Q)ipniQ) belongs to the linear space 

(r-f- — 1 r_ — 1 \ 

Dn + n+ E K + E • 

i+=0 i-=0 J 
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This space is of the dimension + + 1. From the definition of ipn it follows, that 
the functions ipn+i, —rn„ < i < rn+, are the basis of this space. Coefficients Ui{n) of the 
operator Lf are just coefficients of expansion of fipn in the basis of the functions ipn+i- 

The function ipn{Q) can be written explicitly in terms of the Riemann theta- functions. 
For simplicity we present the corresponding formulae for the case r = r+ 



r_. Let a°, 6° be 



a basis of cycles on F with canonical matrix of intersections a° 



W ■ 6° 



0, 



■6° = 5, 



In a standard way it defines the basis of normalized holomorphic differentials ujj{P). The 
matrix B of their 6-periods defines the Jacobian variety J(F) and the Riemann theta-function 
9{z) = 9{z\B). Let us introduce the functions hj, j = 0, . . . ,r — 1, 



where 



e{A{Q) + Zo) nil e{A{Q) + sf) 



9-1 



St = lC-A{Pt)-Y.^(ls). ^ = 0,...,r-l. 



(2.40) 



(2.41) 



i-i 



s=l 

i 



Z, = Z, + Y.A{p-)-Y.A{Pt) Zo = /C-^A(7.), (2.42) 

i=0 i=l s=l 

and A(Q) is a vector with the coordinates Ak{Q) = oJk- Let dVt^^^ be a unique meromorphic 
on F, with simple poles at P^ with \ 

conditions 



differential on F, with simple poles at P^ with the residues =f1, and normalized by the 



h 

The coordinates of the vector of its 6o-periods [/'^"^ equal 



(2.43) 



(2.44) 



Lemma 2.2 The function ipn{Q) equals 

' '^^^(A(Q) + Z,)e(A(P,+) + A;f/(o) + Z,-' ^ ' ' ' ^ 



For the proof of ()2.45p it is enough to show that a function defined by the formula is single- 
valued on F and has all the required analytical properties. 

Corollary 2.1 The coefficients of commuting rank 1 generic operators are quasi-periodic 
meromorphic functions of the variable n. 
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2.4 Rank > 1. The case of separated infinities. 

From the construction of formal common eigenfunctions of the operators L and A it follows, 
that for almost all values of E the operator A{E) is diagonalizable. We call a set of the multi- 
plicities yU = (/ii, . . . , fik) of the eigenvalues of A{E) vector rank of the commuting operators. 
Note, that rank of commuting differential operators with scalar coefficients is always scalar 
(fc = 1), and is a common divisor of orders of the operators |2|. The notion of vector rank in 
the classification problem of commuting differential operators with matrix coefficients was 
introduced in b For commuting operators of a vector rank jj, the characteristic equation 
has the form 

k 

det{w - A{E)) = l[R';^{w,E). (2.46) 

1=1 

We would like to emphasize that conditions ()2.22[ I2.29jl , which are sufficient for simplicity 
of the eigenvalues of the operators A{z),A~{z-) in the formal neighborhood of the infinity 
are incompatible with auxiliary linear problems for the 2D Toda lattice equations. Indeed, 
from the equations 

[L,U] = [A,U]=0, (2.47) 
where L, A h Lj have the form (|1.6jl h (|1.7jl it follows that 

UN+{n + 1) = UN+{n), UN_in)c{n — N^) = UN^{n — l)c(n), 
VM+{n + 1) = VM+{n), VM^{n)c{n - NJ) = VM+{n - l)c{n). 

The later equations imply that vf, given by ()2.7| I2.8jl . satisfy the equations 

vt = v+, v- = v-, (2.48) 

In that case, a priori there are no obstructions for the existence of multiple eigenvalues for 
the operator A{z) or for the operator A^{z^). 

Commuting operators such that the set of eigenvalues of A{z) does not intersect with 
the set of eigenvalues of A~'{zJ) are called operators with separated infinities. Commuting 
operators have separated infinities if at least one of the conditions ()2.38p . i.e. one of the 
conditions (z) m+ra_ ^ m_n+, (ii) 7^ v~ , is satisfied. 

The set of multiplicities /x^ of distinct eigenvalues of the operators A{z), A~{z-) is called 
type of the pair of commuting operators with separated infinities, 

E l^t = (2-49) 

Here I± are finite sets, which parameterize distinct eigenvalues of the operators A{z) and 
A~{z^). Distinct eigenvalues of these operators correspond to distinct "infinity" points on 
the components Fj of the affine spectral curve, defined by the equations Ri{w,E) = 0. A 
simple counting of degree of the divisor of common eigenfunctions (see details below in the 
proof of the first statement of Theorem 2.3 ) shows that there is no component Fj, which is 
compactified by the points corresponding to eigenvalues of only one of the operators A{z) 
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or A (-2-). The multilicity of an eigenvalue of A{E) is constant on each of the components 
Fj. Therefore, 

(i) for each index i± there is at least one index jz^, such that /i^ = /i^. 

As it seems to the authors, there are no more constraints on the types of commuting operators 
with separated infinities, i.e. for any set of positive integers /i^, satisfying equation ()2.49|) 
and the constraint (i), there is a pair of commuting operators of the form ()2.2p with separated 
infinities which has this set as its type. 

Below, we prove the statement for the types of the form (r, r). This is the type of 
commuting operators of the form (|2.2p with separated infinities, which have equal common 
divisors of positive and negative orders, (r = r+ = r_), having the maximum possible scalar 
rank n = r. 

Lemma 2.3 Let a pair of irreducible commuting operators of the form ^2.^) . where (r = 
r+ = r_), have rank r. Then there is a unique Laurent series 

v+{x) = + 0(x-'"++^)) (2.50) 

such that the equations 

= A^{z) = v{z')^{z) (2.51) 
have a solution "^{z) of the form 

= [l + flUn)z'^ . (2.52) 

The space of solutions of equations ^2.51\) in the space L{z) of Laurent series is spanned by 
the series 'i/{e''z),e^ = 1. 



Almost identical statement holds for eigenvalues and eigenvectors of the restriction of the 
operator A on the space L(-2-). 

Let us keep the same notation F for the curve defined by the equation R{w, E) = 0, 
where R{w, E) is a root of degree r from the characteristic polynomial 

det{w-A{E)) = R\w,E). 

From Lemma 2.3 and its analog for the operator A^{zJ) it follows, that at the infinity we 
have 

71+ —1 

n (w-v'[^zk) 



R{w,E) 



k=0 



n_ — 1 



fc=0 



Yl (w -V (Zk-) 



(2.53) 



where the products are taken over all the roots of ^ of the orders n+ h and n_ , respectively. 

To every point of F there corresponds r-dimensional subspace of common eigenvectors 
of the operators L and A. Let us fix a basis ip\Q) in this subspace by the conditions 



iPn{Q) = ^in, i,n = 0, 



,r — 1. 



(2.54) 
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For all n the functions "ipniQ) are meromorphic functions. Their form in the neighborhood 
of the infinity can be found with the help of basic series '$'^{e^z±). The same as in the 
continuous case 0, singularities of the vector- function ipn{Q) = {'^niQ)} in the affine part 
of the spectral curve are described by matrix divisors. 

In the general position, when the poles of ipn are simple, the corresponding matrix divisor 
D = {7s, as} is just a set of distinct points 7^, and a set of r-dimensional vectors = {a*}, 
defined up to a factor as — > Aa^. The points 7^ are poles of ipl^, and the parameters as define 
relations between residues 

aiies^^i^iiQ) = aiies^^Q). (2.55) 

In m El the sets (7, a), were called Tyurin parameters. According to [2^], in a generic 
case they uniquely define a framed stable vector bundle £ over F of rank r and degree 
Ci(det £) = rg. 

Theorem 2.3 Let a pair of irreducible commuting operators with separated infinities sat- 
isfies the conditions and has the rank r = = r^. Then: 

(1) the spectral curve T, defined by the characteristic equation ^2. 5'-^^) is irreducible. It is 
compactified at the infinity by two points , in the neighborhoods of which local coordinates 
are z± = E-i/"±; 

(2) the vector-function ipiQ), whose coordinates are common eigenfunctions 

LC(g)=^^;(g), A^^iiQ)=wi;l,{Q), Q = {w,E)eT, 

normalized by the conditions \2. 54\ ), is meromorphic on F; its matrix divisor outside the 
punctures does not depend on n. In the general position, when the spectral curve is 
smooth, the divisor of poles D = {73,03} is of degree gr, where g is the genus ofV. 

(3a) in the neighborhood of P^ the function iplr+jiQ) ^'^^ ^he form 

^l+i = z-\l + Oiz+)), 

^U, = 0{z-'+'), t>j. (2.56) 
(3b) in the neighborhood of P~ the function i/jl^_^_j{Q) has the form 

= Oizl^'), z<j, 

= Oizt), i>j. (2.51) 

The inverse statement holds. For each generic set of points 7^ and vectors as = (a^), s = 
1, . . . ,gr, i = 0, . . . ,r — 1 there is a unique set of functions ipn{Q) such that: 

(i) ipniQ) outside the punctures has at most simple poles at 7^; their residues at these 
points satisfy relations (I2.55jl : 

(m) TpniQ) at the neighborhoods of the punctures P^ has the form ()2.56l2.57jl . 



preliminary draft 



15 



6 (pespajisi 2008 r. 13:53 



Theorem 2.4 Letip^Q) = {i'niQ)} the functions, defined above by the data {T, P"^, D = 
{ls,C(s}}- Then for each function f G ^(F, P='=) there is a unique difference operator Lf 
(with coefficients independent on Q) such, that 

Lf^\Q) = fiQWiQ). 

If the function f{Q) has poles of orders n+ u n_ at the punctures then the operator Lf 
has the form ^2.^) . where r = r+ = r-. 



Example, r = 2, g = 1. Without loss of generahty we may assume that two punctures on 
the eUiptic curve with the periods {2uj,2u') are points ±zq. Let us choose the vectors ai,a2 
in the form as = (a^, 1). 

From (I2.5eil2.57ll it follows that V'L 

can be represented in the form 



a{z - zo)a{z - 7i - 72 - {2n - l)zo) 
a{z - ^i)(t{z - 72) 



crjz + zq) 
a{z - zo) 



(2.58) 



where a{z) = a{z;2uj,2uj') is the Weierstrass sigma-function. Similar expression for ^/'gn has 
the form 



cr{z - 7i - 2nzQ) (7(2; - 72 - 2nzQ) 



Ciz + Zq) 
a{z - Zo) 



o-(2;-7i) o-(2;-72) 
Conditions ()2.55|1 allow us to express the parameters -B„, Cn in terms of A, 

(t(72 + (2n - l)2o)o-(7i - ^0) 



(2.59) 



Bn = aiAr, 



cririi - -f2)cr{2nzo) 
cr(7i + {2n - l)2;o)cr(72 - 2:0) 



(2.60) 



(2.61) 



^■(72 - 7i)o-(2nzo) 

In the neighborhood of Zq the function ip2n has the form {z — -2o)~". That implies an explicit 
formula for An 



(j(2nzo)o-(7i - 72) 



(ai - a2)fT"(2zo)a((2n - 1)zq + 72)^^(2^ - 1)2:0 + 7i) ' 
In analogous way we find 

,0 ( \ _ A' ^(^ + ^0)^^(2; - 7i - 72 - (2n + l)zo) \(y{z + zo) 

r2n+l[^) — At 



(2.62) 



^2n+l(^) 



a{z - 7i)(t(2; - 72) 
^ D/ (^{z - 7i - 2^2:0) , a{z - 72 - 2nzo) ~ 



(T(Z-7lJ 



a(z-72) 




where 



^/ ^ -1 o-(72 + (2n + l)zo)ff(7i + ^0) 
" " a(7i - 72)a(2nzo) 



(2.63) 
(2.64) 

(2.65) 
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^, = ^(7i + (2n+l).oM72 + .o) ^^ ^^^ 

or(72 - 7i)cr(2n2;o) 

= (2.67) 

w (^(72 + (2n + l)2;o)c^(7i + (2ri - l)^o)t^(7i + 2:0) 

aia{zo - 7i) 

w/ a(7i + (2n + l)zo)o-(72 + (2n - l)^o)a(72 + ^0) „„x 

-/„ = 7 ^ • (2.69) 

02^(^0-72) 

Similar expressions in terms of the Riemann theta-functions can be found in general case. 

Corollary 2.2 Coefficients of the operators Lf, defined by the previous theorem, are quasi- 
periodic functions of the variable n. 

Remark. Note that the functions ?/'2n ^ind ip2n+i iii the continuous limit 

— ^ 0, n 00, nzo x, 

converge to different functions of the continuous variable x. Probably, that explains why 
in the classification problem of commuting differential operators there are no constructions 
without functional parameters. 



2.5 Rank >1. Glued infinities. 

Now let us consider commuting operators of the form p.2|) having maximum possible rank 
/ = r++r_. In this case the formal eigenvalues of the operator A on two infinities L± coincide. 
An eigenvalue of the operator A{z) may coincide with an eigenvalue of the operator A~{z-) 
only when m+ = m_ = m, n+ = n_ = n, i.e. the case of full or partial glued infinities may 
occur only in the classification problem of commuting operators of the form: 

Nr+ Mr+ 

L= <n)T\ A= <n)T\ (iV,M) = l, (2.70) 

i=—Nr- i=—Mr— 

We call a type of commuting pair of operators a set (yU^| /i~), where the index i is in a set 
J, parameterizing all the distinct eigenvalues of A{z) h A~(2;_), and jjif are multiplicities 
of the corresponding eigenvalues. Note, that the type of operators with separated infinities 
can be identified with the particular case of the later definition in which / is a union of sets 
/±, which parameterize indices i± in ()2.49jl . and all the pairs have either the form (/i^| 0) or 

(0|/i7). 

We conjecture that there are no restrictions on the possible types, when the numbers 
/if > are both non-zero. If the type contains a pair of the form (/i| 0) then, as before, 
it should contain also the pair of the form (0|/i). Complete solution of the classification 
problem of commuting difference operator requires a construction of commuting operators 
of all the possible types. We leave the problem open, and consider below a construction of 
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commuting operators of the type that contains only one pair (r_|_| r_). As it was mentioned 
above, that is the case of commuting operators of the maximum possible rank / = r+ + r^. 

Direct spectral problem. From the construction of the formal eigenfunctions at the infin- 
ity it follows, that the necessary conditions for the case of the maximum rank are equations 
vf = v~ = V, i.e. the operators have the from (j2.7()|) . and their leading coefficients in the 
gauge UNr+ = 1 have the from 



U- 



-Nr. = h ^{n — NrJ)h{n), VMr+=v, V-Mr. = vh ^{n — Mr^)h{n). (2-71) 



As above, we keep the same notation V for the curve defined by the equation R{w, E) = 0, 
where R{w, E) is a root of degree I from the characteristic polynomial 

<lei{w-A{E)) = R\w,E). 

At the infinity we have the expansion 



R{w,E) 



■N~l 



Yliw- v{Zk)) 



k=0 



(2.72) 



where the product is taken over all the roots of E~^ of degree A^. Hence, in the case of the 
maximum rank the spectral curve is compactified at the infinity by one smooth point, and 
as a corollary, it is irreducible. 

To every point of T there corresponds /-dimensional space of common eigenfunctions of 
the operators L and A. Let us fix the basis ip^Q) of this space by the conditions 

ijliQ) = Sin, -r.<i,n<r+. (2.73) 

Note, that the choice of the interval for the parameters n, used for the normalization is 
crucial for closed description of analytical properties of the basis functions at the infinity. 

Theorem 2.5 For a generic pair of commuting rank I operators, the common eigenfunc- 
tions ipl^, normalized by ^2. 7,^) . satisfy the following conditions: 

1°. On the affine part of the spectral curve T the functions ipl^ have gl poles 7^ independent 
ofn, at which the equations 

hold. 

2°. In the neighborhood of the "infinity" point Pq the vector-row ipn = {V'n} has the form 

= {T.Un)z'^ *o(n,z), z-^ = E. (2.75) 
Here Cs{n) = '^'^^ vector-rows, 
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\l/o(^5-2) is the Wronsky matrix '^Q\n, z) 
0* of the solutions of the equation 



[z), — r_ < i,j < r^, defined by the basis 



n)(pn+i = z 



(2. 77) 



whose coefficients fi{n) are polynomial functions of the coefficients of the operator L. The 
solutions 0* are normalized by the equation 



*o(0,^) = l. 



(2.78) 



Proof. Let \l/(n, Q), Q G F, be the Wronsky matrix ^^\n,Q) = ipl^_^_j{Q), — r_ < i,j < r+ . 
In the neighborhood of Pq, where z = is a local coordinate, it can be considered as a 

function of the variable z, i.e. \l/(n, 2;). Let us define polynomials 0^ of the variable z~^, by 
fixing their asymptotic behavior as 2; ^ 0. 

Lemma 2.4 In a generic case there are unique functions (j)n{z), that are holomorphic on 
the extended complex z-plane everywhere except at z = 0, and such that in the neighborhood 
of z = the vector-row (f)n{z) = {<Pn{z)) has the form 



r„(z)^(n, 2;), 



(2.79) 



where rn{z) is a vector-row holomorphic in the neighborhood of z = 0, and such, that its 
evaluation at z = equals 

<(0) = (2.80) 

The problem of construction of 0„(2:) is a standard Riemann problem. Let us choose a small 
neighborhood of z = 0. Then we define vector- functions 0„ n r„ as holomorphic outside and 
inside the neighborhood, respectively, and such that on the boundary of the neighborhood 
the equation (12.7911 holds. If the change of the argument of the determinant of \l/(n, z) along 
the contour equals zero, then the Riemann has a unique solution if we fix a value of the 
vector-function at some point. Therefore, for the proof of the Lemma it is enough to show 
that the determinant of \l/(n, z) is holomorphic in the neighborhood of z = and in general 
does not vanish at 2; = 0. The later statement follows from: 



Lemma 2.5 In the neighborhood of the infinity the matrix function 

X{n,Q)=^{n+l,Q)^>~\n,Q) 

has the form 



(2.81) 



( 



X{n, z) 








1 







1 





V X-r_ (n, Z) X-r^+l (n, Z) Xr_+2 ijl, z) 

Xiin, z) = z'^5i,o - fiin,z), 
where fi{n, z) are regular series of the variable z. 








■ Xr+-i{n,z) ) 



(2.82) 



(2.83) 
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Proof. The matrix A" (n, Q) does not depend on the choice of the basis functions Therefore, 
the asymptotic behavior of X{0,Q) in the neighborhood of the infinity, can be found with 
the help of the formal solutions defined in the subsection 2.2. Equations ()2.27|2.;Mjl imply 
that the Wronsky matrix \l/oo(0, z), defined by the formal solutions, has the following block 
form 

where A{z),D{z) are matrices of dimensions (r_ x r_) and (r+ x r+), respectively. All the 
matrices A, B, C, D are regular series in the variable z. The leading term of D{z) is a lower 
triangular matrix with 1 on the diagonal. The leading term of the series A{z) is an upper 
triangular matrix. Hence, the inverse matrix has the form 



' ^ \zCi{z) Di{z) 



(2.85) 



The leading terms of the regular series Ai,Bi,Ci, Di equal 

Ai(0) = A-i(O), 5i(0) = -A-i(0)5(0)D-i(0), 

Di(0) = D-i(O), Ci(0) = -D-i(0)C(0)A-i(0). (2.86) 

The series ipl^ have the form z~^5l + P{z), where /* are regular. Therefore, the last row of 
the matrix A'(O), which is equal to ipr^^iQ), has the form 

X'-+-^^\{]) = z-^5l + r{z). (2.87) 



In the construction of the formal solutions index n can be replaced by n — tiq. That shift 
does not change X{n). Therefore, the last row of the matrix A'(no) has the same structure 
as for A'(O), and the Lemma is proven. 

Normalization (j2.73j) is equivalent to the equation \E'(0,2;) = 1. From (|2.82j) it follows 
that 

n—1 n—1 

det^(n,2) = n detX{n,z) = (-1)" [] f-rAn,z). 

m=0 771=0 

Hence the determinant is holomorphic in the neighborhood of 2; = 0, and in general does not 
vanish at z = 0. That completes the proof of the Lemma 2.4. 

Note, that by their construction the functions 0^ are entire functions of the variable z~^. 
The matrix \l/(n, 2;) is meromorphic in the neighborhood of 2; = 0. Hence, 0^ at 2 = has 
the pole of finite order and, therefore is a polynomial function of z~^. 

Let \E'o(n, 2) be the Wronsky matrix of the functions 0^(2), i.e. \E'-^'*(n, 2) = 0^_,_^(2). 

Lemma 2.6 In the neighborhood of 2 = the matrix function ^'o has the form 

*o(^, z) = R{n, 2)^(n, 2), (2.88) 

where R{n, z) is a matrix function holomorphic at z = 0, and such that R{n, 0) has the block 
form 

i?(n, 0) = ^°^) (8.89) 
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where R+ is a lower-triangular, and R- is an upper triangular {r± x r-t) matrices with 1 on 
the diagonals. 



Proof. From the definition of (^ri it follows that j-th row Rj of the matrix R for j > equals 



Rj{n, z) 



z) \Y '^{^ + h z) 

i=0 



(2.90) 



Equations (l2.8()|01?Hll imply that Rj is regular. Moreover, coordinates Rj{n, 0) of the vector 
Rj{n,0) might be different from zero only for < i < j. At the same time i?j(n, 0) = 1. 
That proves that i?+ is an upper triangular matrix with 1 on the diagonal. The analogous 
statement for i?_ can be proven by similar arguments if one replaces X by the inverse matrix 
which has the form 



X' 



X- 



( X-r- 
1 



V 



X- 



■r_+2 X-r+3 



1 











Xr+-1 1 \ 






(2.91) 



By the definition \E'(0, z) = 1. Entries of \E'o are polynomials in the variable z~^. Therefore, 
p.88|) implies that \E'o(0, 2;) is a constant matrix, equal to i?(0,0). By induction in j, it is 
easy to show that \Efo(0,-2) satisfies the normalization ()2.78ll . i.e. \E'o(0,z) = 1. 

Let us prove that the functions 0^ satisfy an equation of the form ()2.77|1 . Consider the 
matix 

^o(^ + 1, ^)^o ^(^, z) = Xo{n, z) = R{n + 1, z)X{n, z)R{n, z)'^ . (2.92) 

From ()2.83|1 and the structure of -R(n, 0), it follows that the only entry of the last row in 
A'o, having pole at 2; = is A coefficient at the singular term z~^ of its expansion 

equals 1. The matrix is regular for all 2; 7^ 0. Therefore, it is a polynomial function of 
z~^. Hence the last row of Xq has the desired form. 

In order to complete the proof of the second statement of the Theorem it is enough to 
invert equation (j2.88|] . We have 

^ = i?^i^o- (2.93) 

This equation for the row with the index j = gives ()2.75p . where the first factor is the 
Tailor expansion of the corresponding row of the matrix R~^{n, z). 

Using the asymptotic of ipn at the infinity one can find the degree of the matrix divisor 
of its poles using arguments identical to that in the continuous case [2]. 

The Theorem 2.3. assigns for the generic pair of commuting operator of the form ()2.7n|) . 
having rank / = (r+ + rJ) : a smooth algebraic curve F with a puncture Pq; a set of Tyurin 
parameters and a set of / — 1 functions fiiji) of the discrete variable n 



L,A^{F,(7,a),/,°(n)}. 



(2.94) 
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Inverse spectral problem. Let F be a smooth algebraic curve with fixed local coordinate 
z at the neighborhood of puncture Pq- Let us fix a set of function fi^{n), r_ < i < r+. They 
define the Wronsky matrix "^q, \I/o * = 0^+^, < j < r^, corresponding to the solutions of 
the degree Z = r+ + r_ difference equation 

£ /,°(n)0„+, = ;^-Vn, /°^ = 1, (2.95) 

i=r— 

normalized by the conditions 

<Pl = 6i r.<t<r+. (2.96) 



Theorem 2.6 fj^ij/j For a generic set of Tyurin parameters of degree Ig and rank I, i.e. for 
a generic set oflg points 7^ and a set of projective l-dimensional vectors as = (««), r_ < z < 
r+ there is a unique vector-function ipn{Q) with coordinates, which outside the puncture Pq 
have at most simple poles at 7^. Their residues satisfy relations \2.74 ). In the neighborhood 
of Pq the vector-row tpn has the form 

= (^Un)z'^ ^o(ri,^), a = ^0- (^-97) 

Let f G ^(r, Pq) he a meromorphic function on T with the only pole at Pq of order N . Then, 
there is a unique operator Lj of the form 

Nr+ 

Lf= Uiin)T\ UNr+ = l, (2.98) 

i=-Nr- 

such that 

Lf^iQ)=fiQ)m)- (2.99) 

A proof of the theorem is standard. Equation (I2.97jl means that ^ is a solution of the Riemann 
problem on F, in which \E'o is the transition function in the neighborhood of the puncture. 
In the general position the existence and uniqueness of the solution of this problem follows 
from the results of |27[ I28| . or simply from the Riemann- Roch theorem for vector bundles 
(see details in [2]). The same results imply the second statement of the theorem. 



Theorem 2.7 A commutative ring A of operators of the form ^2. 9^) of maximum rank I is 
isomorphic to the ring A(T, Pq) of meromorphic functions on an algebraic curve T with the 
only pole at the puncture Pq. In the general case the isomorphism ^(F, Pq) = A is defined 
by equation ^2. 99\) . where the Baker- Akhiezer function corresponds to a set of the Tyurin 
parameters (7,0). 
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2.6 Discrete dynamics of the Tyrin parameters 



In this section discrete equations for the Tyurin parameters corresponding to commuting 
operators of maximum rank / are derived. According to Theorem 2.6, a smooth algebraic 
curve r with a puncture Pq, a generic set of the Tyurin parameters (7,0) and coefficients 
fiiji) of difference equation ()2.95p define the vector- function ^n(Q)- The corresponding 
Wronsky matrix is denoted by '^{n,Q). The matrix function X{n,Q)^ given by 1)2.811) . has 
the asymptotic (I2IH21 EES!) • Let 

/^W = /i(n,0) (2.100) 

be the evaluations of the regular series fi{n,z) in ()2.83|) at the puncture z = 0. These 
functions of the discrete variable n can be expressed in terms of the initial parameters f^{n), 
and the first coefficients iiiji) of expansion (|2.97|] for ipn- The corresponding formulae are far 
from being effective, because the coefficients iiiji) depend on the initial data (7^, a^, fiiji) 
through the Riemann problem. At the same time, as it will be shown below, there is no 
need for explicit formulae for fi{n). They can be chosen as a part of independent parameters 
defining commuting operators. 

For n 7^ let us denote zeros of det \E'(n, Q) by 7s (n). In general position, when they are 
simple their number equals gl. Let as{n) be the corresponding left zero- vector 

asH*(n,7,(n)) =0. (2.101) 

For n = 0, we set 

7.(0)= 7s, «s(0) = as. (2.102) 
From the definition (12.8111 it follows: 

Lemma 2.7 The matrix function X{n, Q) has simple poles at the points 7s(n) and satisfies 
the relations 

aiin) Tes^^^n)X""\n, Q) = ai{n) res-,4„)A^'"'^(n, Q). (2.103) 
The points 7s (n + 1) are zeros of X{n, Q), 

detA'(n,7s(n + 1)) = 0. (2.104) 

The vector as{n + 1) is the left eigenvector of the matrix A'(n,7s(n + 1)), 

as(n+l);f(n,7s(n + l)) = 0. (2.105) 

Simple dimension counting implies the following statement. 

Lemma 2.8 Let T be a smooth algebraic curve with fixed local coordinate c z{Q) in the 
neighborhood of a puncture Pq. Then for a generic set of data (7s(^), «s(^), fii"^)) there 
is a unique meromorphic matrix function X{n,Q), Q gT, with at most simple poles at the 
pointsPo,'ys, and such that: 

{i) the Laurent expansion of X{n,Q) at Pq has the form ^2.8^^2. 8,^) . where the regular 
series fi{n,z) satisfy ^2. 101 A): 

{ii) the residues of X{n,Q) at 7s satisfy the relations ^2.10,^) . 
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Equations (12 ■1041 I2■lfl5^1 can be regarded as equations that define {'ys{n + + 1)) 

through the matrix X{n,Q). The later is uniquely defined by {'ys{n),as{n), fi{n)). Hence: 

Corollary 2.3 The functions fi{n) and the Tyurin parameters (7, a), which define the ini- 
tial condition ^2.10^) of the corresponding discrete dynamical system, are the full set of 
data, which parameterize commuting operators with the given spectral curve. 



Example g = 1, 1 = 2. Let us consider a pair of commuting operators 

2 3 
L=Y1 Ui{n)T\ A=Y. v,{n)T (2.106) 

i=-2 i=-3 

of the maximum rank / = 2. Their spectral curve T is an elliptic curve. Let 2uj, 2uj' be periods 
of r. Without loss of generality we identify the puncture Pq with the point z = in the 
fundamental domain of z-plane. 

The operators (12.10611 are uniquely defined by Tyurin parameters, and the parameters 
fi{n), i = —1,0, which are denoted below by 

/_i = c„+i, /o = Vn+l. (2.107) 

Our next goal is to find explicit formulae for the coefficients of commuting operators ()2.106j] 
with the help of equations of the discrete dynamics of Tyurin parameters. For brevity, let 
us introduce the notations 

7n=7lH, ln = l2{n). (2.108) 

For the vectors a*, i = —1,0, which are defined up to a factor, we fix the normalization 
under which as{n) are two-dimensional row- vectors with the coordinates 

ai{n) = {all), a2{n) = (all). (2.109) 

According to Lemma 2.8, the data 7^'^, a^'^, c„+i, Vn+i uniquely define the matrix A"^* = 
Xl^Ujz), i,j = —1,0. Let us find its explicit form in terms of the standard Weierstrass 
functions. By definition, Xn has the form 

An elliptic function Xn(^) has poles at 7^'^ and equals — c^+i at 2; = 0. Therefore, it can be 
written as 

xl = -c„+i + (c(^ - 7^) + C(7;^)) + B, (c(^ - il) + Ciin)) , (2.111) 

where ({z) is the Weierstrass C-function. 

The function xi in the neighborhood of the puncture z = has the form = ~ 
Vn+i + 0{z), i.e. 

Xl = -Vn+l + C(^) + ^2 (C(^ - 7n) + C(7n)) + (C(^ " 7n) + C(7')) • (2-112) 
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The sum of residues of an elliptic function equals zero. Hence, 

Ai + 5i = 0, ^2 + 52 = -!. (2.113) 

At the same time 

Ai = alA2, Bi = alB2. (2.114) 
From (12.11 ,312.1 14jl it follows that 



X\ = -Cn+i + -T3t^(C(^-7;i,)-C(^-7D + C(7^)-C(7^)), (2.115) 

xl = C(^)-^n+i + ^^(C(^-7^) + C(7;i)) + 

+^^(C(^-7D+C(7D) . (2.116) 

According to Lemma 2.7, the points Jn+i defined by the equation 

detXn{Yn+i) = x\Yn+i) = 0. Hencc, 



Cn+i = (C(7:+i - 7^) - C(7:+i - 7^) + C(7^) - C(7^)) . (2.117) 

Note, that the sum 7^ + 7,^ = 2c does not depend on n, because 7^ are poles and 'Jn+i 
zeros of an elliptic function. Therefore, 

7n = 7n + c, 7^ = -7„ + c, c = const. (2.118) 

Using p. lisp , equations p.ll5|2.116p can be rewritten as 

Xn = 1 " " 2 - 7n - C) - C{Z + 7„ - c) - C(7n+1 ' In) - C{ln+1 + In)] (2.119) 



a. 



2 



= -Vn+1+C{z)+ ^ " ^ (C(^-7n-c)+C(7n + c)) + 



+ (C(^ + 7n - c) - C(7n - c)) . (2.120) 

If the functions f„ and 7^ of the discrete variable n are chosen as independent parameters, 
then equation p.ll7p 

Cn+l = (C(7n+1 - 7n) " C(7n+1 + In) + ({in + c) + C(7n " c)) (2.121) 

can be regarded as the definition of the variables c„+i. 
From (12.10511 it follows that 

a,{n+l) = -xl{ls{n + l)). (2.122) 
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Equation (|2.12fljl . implies recurrent relations defining a. 



1,2 
n+1- 



a„+i = t;„+i-C(7n+i + c)-^-^(C(7n+i-7n)+C(7n + c)) 



al 



-(C(7n+l+7n)-C(7n-c)), (2.123) 



at 



a: 



n+l = Vrr+l + C(7n+1 " c) + " (C(7n+1 + In) - C(7« + c)) 



al 



(C(7n+l-7n)+C(7n-c)). (2.124) 



As shown above, the arbitrary functions 7^, f „ and a certain constant c define the matrix 
function Xn, and as a corollary, define the coefficients of commuting rank 2 operators, cor- 
responding to an elliptic curve. Each operator corresponds to an elliptic functions with the 
pole at the puncture z = 0. The simplest operator L4 is of order 4 and corresponds to the 
Weierstrass function p{z) = z~'^ + 0(2;^). 

Coefficients of this operator are coefficients of the expansion of ip{z)il)n^^ in the basis 
?/'ri+i\I/~^, 2 <i<2. The expansion is determined by singular parts of the Laurent series at 
z = 0. Let tpm be polynomials in k = z~^, such that ipm'^n^ = "ipm + 0{k~^). Then 



= {-Cn+l,k - Vn+l), ^n+l = (0, 1), = (1, 0). 

The polynomial ipn-i can be found with the help of the equation = X^^'ifn, 
Let be coefficients of the expansions 

Xl = k- Vn+l + C^Z + ■■■ 

Then 

In a similar way ipn-2 can be found. After a long but straightforward calculation we 

L, = Ll- + T + Cn {ail + + , 



2.125) 
2.126) 

2.127) 
2.128) 

2.129) 

find 

2.130) 



where L2 is the discrete Schrodinger operator with the coefficients 

L2=T + Vn + CnT~\ (2.131) 

and the function m„ is given by the formula 

(2.132) 
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Symmetric case. Let the constant c in (j2.118|l be zero, c = 0. Then xli even function 
of z. Hence, = 0. Equation (12.1191) implies 



^12 _ «n«n P'iln) 



e = r \, ■ (2-133) 



Using p.l2H] . we obtain 



e = T7 ^ , f,^^""^ , , = P(7n) - P(7n+l) . (2.134) 

C(7n+1 + In) - C(7n+1 " 7n) - 2C(7n) 

(The last equation can be obtained from the addition formulae for the Weierstrass ^-function, 
or by direct comparison of the poles and the residues of both sides of the equation.) From 
(12.1 2011 it follows that 

C = Piln) ■ (2.135) 

Substitution of the last two formulae into ()2.132p implies that the operator L4 in the sym- 
metric case equals 

L4 = L2 - p(7„) - p(7n-i)- (2.136) 

In the symmetric case the formulae for the coefficients of L2 can by simplified. Let F{u, v) 
be the elliptic function 

F{u, v) = ({u + v)- ({u -v)- 2av) = f ^"^^ . (2.137) 

Then, the formulae (j2.121ll2.124|] for the symmetric case c = can be represented in the 
form 

" - i^(7n+l,7n), (2.138) 



al^^ = vn+i + \ ( F(7„, 7„+i) + ^i±4F(7„+i, 7„) ] , (2.139) 



a, 



2 

n+l 



n+1 - \ fi^(7n,7n+l) - 4^^(7n+l , 7n)) • (2.140) 



The last two equations are equivalent to 

a^+i-a^+i = F(7„,7„+i), (2.141) 

a^+i + a^^i = 2t;„+i + 4^i^(7n+i,7n). (2.142) 



Let us define Sn by the formula 



Then 



(2.143) 



+ = -^ni^(7n-l,7n), = " l) i^(7n-l, 7n) (2.144) 

preliminary draft 27 6 (pespajia 2008 r. 13:53 



and equations (I2.138[ 12.1421) can be rewritten as 

4Cn+l = (4 - l)^(7n+l,7n)i^(7n-l,7n), 



(2.145) 



'5?i-^(7n+l) 7n) ni In+l) ■ 



(2.146) 



The last equation shows, that in the symmetric case the variables 7^, Sn can be chosen as 
independent parameters. Then, equations (I2.131f2.i;ffil2.145|2.146l) give explicit expressions 
for the coefficients of the operator L4, which were presented in the introduction. In the 
similar way one can find coefficients of the operator Aq . 

3 High rank solutions of the 2D Toda lattice equations. 

A key element of the algebro-geometric integration theory of non-linear equations is a con- 
struction of the multiparametric Baker- Akhiezer functions. These functions are defined by 
their analytic properties on a corresponding algebraic curve. Below we define the multi- 
parametric Baker- Akhiezer vector-functions. They are deformations of the eigenfunctions 
of commuting operators of an arbitrary rank. The definitions are different in the cases of 
operators with separated or glued infinities. 



3.1 Separated infinities. 

As it was shown above, operators with separated infinities are defined by their algebro- 
geometric data: a spectral curve with two punctures and a set of Tyurin parameters. There no 
functional parameters in their construction. Functional parameters do appear in construction 
of the corresponding solutions of the 2D Toda lattice equations. It is necessary to stress that 
these functional parameters are functions of continuous variables, but not a discrete one, as 
in the construction of commuting operators with glued infinities. The functional parameters 
define germ functions \E'±(t^,z), depending on the corresponding half of the times of the 
hierarchy of the 2D Toda, and are entire functions of the variable z~^. 

Consider functions '^±{z) of the variable z~^, such that 



Lemma 3.1 Let T be a smooth genus g algebraic curve with fixed local coordinates z± in 
the neighborhoods of two punctures P±. Then for a generic set of the Tyurin parameters of 
degree rg and rank r there is a unique vector-function i'niQ) such that: 

(z) its coordinates ipl^, i = 0, ... ,r — 1, outside the punctures have at most simple 
poles at the points js, where the relations ^2.55\} hold; 

(a) in the neighborhoods of the punctures the function ipn has the form 




(3.1) 



ipkr+j = zl''R±{kr+j,z±)'i'±{z±) , 



(3.2) 



preliminary draft 



28 



6 (pespajisi 2008 r. 13:53 



where R±{n^z) are holomorphic in the neighborhood of z = row-vectors with values at 
z = satisfying the normalization conditions 



Ri{kr + J, 0) = 0, i<j, Ri{kr+j,0) 



« > J, 
i =j. 



(3.3) 



Equations ()3.2|l are equivalent to the Riemann factorization problem on F. The dimension of 
the space of its meromorphic solutions with poles at 7^ can be found from the Riemann- Roch 
theorem. This dimension equals the number of linear equations (12 .551) . 

If the functions ^± = '^±{t^,z) depend on some independent variables = {tf), then 
the corresponding vector-functions ipn are functions of the full set of the variables ipn = 
ipnit^ 1 1' ,Q)- Let us define the dependence of \l/± on t^ such, that the corresponding Baker- 
Akhiezer functions give solutions of the 2D Toda lattice. 

If we restrict ourselves to the construction of the 2D Toda equations (jl.ip . only, then the 
dependence on tf = ^,t{ = rj is defined by ordinary differential equations. Their coefficients 
are the functional parameters mentioned above. 

Let ai(tf),yi{t]^) be a set of arbitrary functions. Then the germ functions \I'± are defined 
with the help of the equations 



^±(0,^) = 1, 



(3.4) 



where the matrices M^^{tf, z) have the form 



0,1 



/ 




1 


■ 


■ 






( 
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■ 


boz- 


1 \ 







ai 


1 ■ 


■ 






hi 
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■ 


























■ 
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■ 1 
















V 


z-^ 
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■ a^^i 


) 




I 


■ 


■ K-i 





/ 



(3.5) 



where 



(3.6) 



Theorem 3.1 The Baker- Akhiezer vector-function tpn, corresponding to a generic set of 
data {T, P±, z±, {'j,a),ai, yi} satisfies the equations 

with the coefficients 
where 

ifkr+i = Vrit^) + In Rt{kr + i, 0, tf, t^). (3.8) 

For each function f G AiV., P^) with the poles of order n± at the punctures there is a 
unique difference operator Lf of the form ^2. ^) with coefficients depending on tf such that 

Lfr = fr- 
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The compatibility condition of equations ()3.7jl is equivalent to the 2D Toda lattice equation 

(HJ. 

Corollary 3.1 The functions ^Pn, given by ( t.V. 7j) . where the second term is defined by the 
evaluation at zero of the regular factor R- in (E^, are solutions of the 2D Toda lattice 
equations. 

Proof of the Theorem is standard and follows from comparison of analytical properties on 
r of the functions defined by the left and the right hand sides of ()3.7I) . The proof does not 
relay essentially on the specific form ()3.5|l of the matrices M^'^. The Theorem remains true 
if the matrices m£'^ are replaced by matrices of the form 

Ml'' = Ml'' + m±{tf), (3.9) 

where entries m± do not depend on z and satisfy the constraints 

m+ = 0, i < j, rnl = 0,i > j. (3.10) 

However, this extension of a set of the germ-functions provides no new solutions of the 2D 
Toda lattice equations. Indeed, expansion ()3.2p and constraints ()3.3p are invariant under the 
transformations 

^± = ^?±^±, R± = R±gl\ (3.11) 

where g^it^) is an upper-triangular matrix, and g+{ti) is a lower-triangular matrix with 1 on 
the diagonal. Therefore, these transformations do not change the Baker- Akhiezer function. 
They correspond to the gauge transformations 

M^^ gl^d,±g^ - gl^Ml'^g^. (3.12) 

In each gauge equivalence class there is a unique matrices satisfying the constraints m± = 0. 

In order to obtain solutions of the full hierarchy of the 2D Toda lattice equations it is 
enough to define a dependence of the germ functions 'if± on the times tp with the help of 
the differential equations 

9,±^± = M°'P^±, ^±(0,z) = l, (3.13) 

where the matrices M^'^{t^,z) are polynomial functions of the variable z~^. Compatibility 
conditions of ()3.13|) for each half of the times t^ or t~ are gauge equivalent to one of the 
r-reducation of the KP hierarchy. Because our mein goal is a construction of solutions of 
equation (II. Ill , we will present an explicit description of the matrices M±^, p > 1, and 
detailed analysis of the corresponding auxiliary soliton system elsewhere. 

3.2 One-puncture case 

In in the previous case of the separated infinities, the dependence of one-puncture Baker- 
Alkhiezer function on the variables t^ is completely determined by the dependence on these 
variables of the germ function \l/o(^, t, z). The later is defined with the help of the equations 

9,±^o(^, t, z) = M^P{n, t, z) ^o{n, t, z) , (3.14) 
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where the matrices M^^{n,t, z) are polynomial functions of z ^. Unlike the case of the 
separated infinities, even for the construction of solutions to the 2D Toda lattice equation 



the matrices should satisfy the compatibility condition for equations ()3.14|] for p 
and the difference equation 



"^oin + 1, t, z) = Xo{n, t, z)'^o{n, z) 



(3.15) 



which is equivalent to the definition of \l/o as the Wronsky matrix, corresponding to solutions 



of ESS). The matrix A'n 



An 



/ 








1 







1 



1, has the form 








X- 



r_+2 



(3.16) 



where 



rO,l 



B{n,t)XQ 



-1 



(3.17) 
(3.18) 



X° = ^-^5,o-/.°(n,t). 

Let M^'^ have the form 

M^^ = Xo + A{n,t), M'l 

where A and B are diagonal matrices 

A = diag{a{n—r^, t), . . . , a(n+r+ — 1, t)}, B = diag{b{n—r_, t), . . . , b{n+r^ — l,t)} (3.19) 

Then, compatibility of equations (|3.14jl for p = 1 and (|3.15jl is equivalent to compatibility 
of the system 

= (j)n+i + a{n, t)(f)n, dt~<Pn = b{n, t)0„_i, 

and equation (j2.95j) . Therefore, the functions y„(t) such that b{n,t) = e^"~^""% are the 
solution of some reduction of the 2D Toda lattice equations onto stationary points of a 
linear combination of the hierarchy flows, corresponding to the times . . . , t^. 

Let us fix a solution ?/„ of this reduction and denote the corresponding solution of the 
linear system ()3.14[ I3.15|) by '^o{n,ti ,ti , z). Then the following statement is valid. 

Theorem 3.2 (\21^) For each smooth genus g algebraic curve T with fixed local coordinate z 
in the neighborhood of a puncture Pq, and each generic set of the Tyurin parameters (7, a) 
of degree Ig and rank I, there is a unique vector function ipnitf ,ti ,Q) , with coordinates 
which outside the puncture Pq have at most simple poles at the points 7^. Their residues at 
these points satisfy relations \2.74 )- In the neighborhood of Pq the row-vector ipn has the 
form 



J2 Un, tt, ti)z' ^o{n, tt, t^,z), a = S'o- 

KS=0 ) 



The functions ipn satisfy the equations 



d^+lpn = + {dt+^n) V^n, d^~1pn = (c'^" ^„-l, 



(3.20) 



(3.21) 
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where (fn is given by the formula 

y^n = yn{tt,t^)+in {l + ^[-'\n,tt,t^)), (3.22) 

in which is the coordinate with the index i = —1 of the vector in ^S.20\) . 

Example. In the case of rank 1 = 2, r± = 1 the germ function \I/o is defined by any solution 
of the one-dimensional Toda lattice equation 

= ey^-y--^ - ey^+'-y-, (3.23) 

and has the form 

^0 = $(n,t,2)e''^"\ x = tt + t^, t = tt-t^, (3.24) 
where $ is the Wronsky matrix of the auxiliary linear problem for (|3.23|1 . 

3.3 Deformations of the Tyurin parameters 

The reconstruction problem of the Baker- Akhiezer vector-function is equivalent to the linear 
Riemann problem. As it was mentioned above, this problem in generic case can not be solved 
explicitly. At the same time the corresponding solutions of the 2D Toda lattice can be found 
with the help of the equations for the deformation of the Tyurin parameters. 

Let \l/(n, t, Q) be the Wronsky matrix whose rows are the vector Baker-Akhiezer functions 
ipn+j(t,Q)- The corresponding deformation of the Tyurin parameters is defined as follows. 
In the generic case det \l/(n, t, Q) has gl simple zeros 7s(n, t), and as{n,t) are defined as the 
corresponding left null-vector 

a,(n,t)^(n,t,7,(n,t)) = 0. (3.25) 

Difference equations describing the dynamics of Tyurin parameters with respect to the vari- 
able n, were obtained in Section 2.5. The equations for the continuous deformations follow 
from the earlier results of the authors |2T] . 

Let us consider the logariphmic derivative of \l/ 

d^±^ = M^^. (3.26) 

It is a meromorphic function on T and outside the puncture has simple poles at the points 
7s = 1s{n,t). In the neighborhood of 7^ it has the form 

M= ' +^is + 0{z-z{^,), (3.27) 

where is a vector-column. (For brevity we skip indices p, ±.) The first two coefficients of 
expansion 1)3.27^ define the equations of deformation with respect to the variable t = t% 

dt z{'^s) = -Tr (m^Os) = -("s^is), dta^ = -a^fis + ^.s^s- (3.28) 
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Here Kg is a constant. Equation (|3.28jl is a well-defined dynamical system on the space of 
Tyurin parameters, which is the symmetric power S^^ (V x CP^~^^ . 

The compatibility condition 

dtXn = Mn+lXn - XnMn (3.29) 

of the linear problems 

= Xn-^n, dt-^n = M„^„, (3.30) 

is equivalent to a well-defined system of non-linear equations on singular parts of Xn h M„ in 
the neighborhood of the puncture. Here and below \l/„ = \l/(n, t,(5), A'„ = X{n,t,Q), Mn = 
M{n,t,Q).) 

Discrete analog of the Krichever-Novikov equation. Let us consider as an instructive 
example the nonlinear equations in the case / = 2 and g = I- Recall, that in this case the 
coefficients of the linear system defining the germ function $ in (I3.24jl have the form 



^°=( ,0°'^ ,0 )^ Mn=(~'^,f- 1), k = z-\ (3.31) 



The Lax equations for this systems are equivalent to the equations of the one-dimensional 
Toda lattice. 

The leading terms of the Laurent expansion of the "dressed"matrices Xn (see. 1)2.110^ ) 
and Mn have the same form ^3.31^ but with different functions c„,f„. In particular, the 
matrix M„ in the neighborhood of z = has the form 



M =f 2.„-A:, 2 
'^'^ I -2c„+i, k 



j + rrink^^ + 0{k-^), k = z'K (3.32) 



Equations ()3.29|) imply 

Cn+l = 2Cn+liVn+l - Vn), Vn+1 = 2(c„+2 " C„+i) + - mf^^. (3.33) 

Additional terms can be expressed through c„, Vn and the Tyurin parameters 7^, a^. Our 
goal is to get a closed system of equations using the dynamics of Tyurin parameters. 

For simplicity we consider the symmetric case, in which the constant c in p.llSp equals 
zero, c = 0. From the definition of M^, it follows that 

Mf = -c„+i + A'2i, Mf = + A'f . (3.34) 

Hence, 

<' = e'' = p(7n). (3.35) 
Substitution of this formula in b ()3.33|) gives 

Vn+l = 2(c„+2 - Cn+l) + p{ln) - p{ln+l) ■ (3.36) 
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Equation ()3.28jl implies 

a + CL 



% = -res^„M„ = 1 . 3.37 



Therefore, % can be identified with the variables defined in ()2.143|) . and equations ()2.1451 
12.14611 can be represented in the form 

Acn+i = (7^-l)F(7„+i,7„)F(7„_i,70, (3.38) 

2Vn+l = 7n-^(7n+l5 7n) - 7n+l-^(7n,7n+l)- (3.39) 

Let us present two identities, which will be used below. 

du\nF{u,v) = -F{v,u), (3.40) 

dJnF{u,v) = -F{u,v) + 2C{2v) -ACiv), (3.41) 

where the elliptic function F{u,v) is given by ()2.137p . The identities can be verified by 
comparing of singularities of their right and left hand sides. Substitution of ()3.38I3.39|) into 
the first equation (I3.33jl . gives with the help of (|3.4()[ imT|l the equation 

7„ = (7. - 1) (^(7n, 7n+l) + ^(7n, 7n+l)) , (3.42) 

where 

V{u, v) = Ciu + v) + Ciu -v)- C(2m). (3.43) 

In the same way it is possible to check that the substitution of (j3.38l3.39jl into (j3.36ll gives 
the same system (I3.42jl . 

Equations (I3.42jl are Hamiltonian system of equation with the Hamiltonian 

H = J2^n (sh-2 (Pn/2)) + In (p(x„ - Xn-i) - p(a;„ + a;„_i)) . (3.44) 



This system was obtained by one of the authors in fSOj, as the solution of reconstruction 
problem of an integrable system corresponding to a given family of the spectral curves. 
Inverse problems of this kind are natural in the framework of the Seiberg-Witten theory in 
which families of curves parameterize moduli of nonequivalent physical vacuum states in 
supersymmetric gauge models. 

In [nOl the system ()3.42|1 was called elliptic analog of the Toda lattice. After a suitable 
change of variables it coincides with one of systems obtained in IHl] in the framework of 
solution of the classification problem of integrable chains. In (20] the system ()3.42|1 was 
identified with the pole system, describing solutions of the 2D Toda lattice that are elliptic 
in X. An appearance of the same system in the theory of rank 2 solutions of the 2D Toda 
lattice equation came for authors as a complete surprise. 
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